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Abstract 

We review the properties and applications of binary and millisecond pulsars. Our knowl- 
edge of these exciting objects has greatly increased in recent years, mainly due to successful 
surveys which have brought the known pulsar population to over 1300. There are now 56 
binary and millisecond pulsars in the Galactic disk and a further 47 in globular clusters. 
This review is concerned primarily with the results and spin-offs from these surveys which 
are of particular interest to the relativity community. 
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1 Preamble 



In the 34 years that have elapsed since the discovery [ 115 | of pulsars, rapidly rotating 
highly magnetised neutron stars, the study of these fascinating objects has resulted in 
many applications in physics and astronomy. Striking examples include the confirmation 
of the existence of gravitational radiation as predicted by general relativity ||244| , P45|| 
and the first detection of an extra-solar planetary system [|274| , ^f. 




Figure 1: The numbers and locations of the various types of radio pulsars known as of 
December 2000. The large and small Magellanic clouds are denoted by LMC and SMC. 

The diverse zoo of radio pulsars currently known is summarized graphically by the 
Venn diagram in Fig. 0. Many new binary systems containing neutron stars are being 
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discovered as a result of the latest generation of pulsar surveys. This review is concerned 
primarily with the results and spin-offs from these surveys which are of particular interest 
to the relativity community. 



What's new in this review? 



Since the first version of this article was written back in 1997/8 ||146|| a number of pulsar 



surveys using the Parkes radio telescope ||173|| have discovered almost 700 pulsars. As a 



result, the sample size is now double what it was in 1997. Many of the exciting new discov- 
eries from these searches are discussed in this review. Up-to-date tables of parameters of 
binary and millisecond pulsars are included as an appendix. Several new sections/figures 
have been added and existing sections reworked and modularized to make the review 
more self-contained and (hopefully!) easier to read in an html setting. We begin in §0 
with an overview of the pulsar phenomenon, the key observed population properties, the 
origin and evolution of pulsars and an introduction to pulsar search techniques. In §|^, 
we review present understanding in pulsar demography, discussing selection effects and 
the techniques used to correct for them in the observed sample. This leads to robust 
estimates of the total number of normal and millisecond pulsars (§ p.3|) and relativistic 
binaries (§ |3.4| ) in the Galaxy and has implications for the detection of gravitational ra- 
diation from these systems. We discuss pulsar timing in §H. One apphcation of these 
exceptional clocks, a sensitive detector of long-period gravitational waves, is discussed in 
§||. We conclude with a brief outlook to the future in §|[ 



2 An Introduction to Pulsar Astronomy 

Many of the basic observational facts about radio pulsars were established shortly after 
their discovery ||115|| by Bell and Hewish in 1967. In the intervening years, theoretical 
and observational progress has flourished. Although there are many remaining questions, 
particularly about the emission mechanism, the basic model has long been established 
beyond all reasonable doubt viz: Pulsars are rapidly rotating, highly magnetised neutron 
stars formed during the supernova explosions of massive (~ 5-10 Mq) stars. In the 
following, we discuss the basic observational properties most relevant to this review. 

2.1 The Ughthouse model 

An animation showing the rotating neutron star or "lighthouse model" of the basic pulsar 
phenomenon is shown in Fig. |^. As the neutron star spins, charged particles are accel- 
erated out along magnetic field lines in the magnetosphere (depicted by the light blue 
cones). This acceleration causes the particles to emit electromagnetic radiation, most 
readily detected at radio frequencies as a sequence of observed pulses produced as the 
magnetic axis (and hence the radiation beam) crosses the observer's line of sight each 
rotation. The repetition period of the pulses is therefore simply the rotation period of the 
neutron star. The moving "tracker ball" on the pulse profile in the animation shows the 
relationship between observed intensity and rotational phase of the neutron star. 

Neutron stars are extremely stable rotators. They are essentially large celestial fly- 
wheels with moments of inertia ~ lO''^ g cm^. The rotating neutron star model, indepen- 
dently developed by Pacini and Gold in 1968 ||194| , |103|| , predicts a gradual increase in the 



pulse period as the outgoing radiation carries away rotational kinetic energy. This model 
became universally accepted when a period increase of 36.5 ns per day was measured for 
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Figure 2: The rotating neutron star (or "lighthouse") model for pulsar emission. Click 
here to see the movie in action. Animation designed by Michael Kramer. 



the pulsar in the Crab nebula ||2U9|| , enabling Gold | ]104| to show that a rotating neutron 
star with a large magnetic field must be the dominant energy supply for the nebula. 



2.2 Pulse profiles 

Pulsars are weak radio sources. Measured flux densities, usually quoted in the literature 
for a radio frequency of 400 MHz, vary between 0.1 and 5000 mJy (1 Jy = 10"^^ W m~^ 
Hz~^). This means that, even with a large radio telescope, the coherent addition of many 
thousands of pulses is required in order to produce an integrated proflle. Remarkably, 
although the individual pulses vary dramatically from pulse to pulse, at any particular 
observing frequency the integrated proflle is very stable. The pulse proflle can thus be 
thought of as a flngerprint of the emission beam. The animation in Fig. ^ shows a 



Figure 3: Single pulses from PSR B0329-h54. Click here to see the movie in action. 

sequence of consecutive single pulses from PSR B0329+54j^ one of the brightest pulsars. 
This pulsar is seen in the animation to stabilise into its characteristic 3-component form 
after the summation of a number of seemingly erratic single pulses. Stabilization time- 
scales are typically several hundred pulses ||112|| . This property is of key importance in 



pulsar timing measurements discussed in detail in 

^Pulsars are named with a PSR prefix followed by a "B" or a "J" and their celestial coordinates. Those 
pulsars discovered prior to 1990 are usually referred to by their "B" names (Besselian 1950 system). More 
recent discoveries are only referred to by their "J" names (Julian 2000 system). 
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Fig. ^ shows the rich diversity in morphology from simple single-component profiles 
to examples in which emission is observed over the entire pulse. The astute reader will 
notice two examples of "interpulses" — a secondary pulse separated by about 180 de- 
grees from the main pulse. The most natural interpretation for this phenomenon is that 
the two pulses originate from opposite magnetic poles of the neutron star (see however 
|174| |). Since this is an unlikely viewing angle we would expect interpulses to be a rare 



phenomenon. Indeed this is the case: the fraction of known pulsars in which interpulses 
are observed in their pulse profiles is only a few percent. 




(d) B0525 + 21 (e) J1012 + 5307 (f) B1B31-04 




(g) B1937 + 21 (h) J2124-3358 (i) J2145-0750 

j 
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Figure 4: A variety of integrated pulse proRles taken from the available literature. Refer- 
ences: (a,b,d,f: / [105|/ ); (c: J3^/); (e,g,i: !\139{] ); (h: ^^). Each profile represents 360 degrees 
of rotational phase. These profiles are part of a database of over 2600 multi- frequency 
pulse profiles for over 600 pulsars that is available on-line 




(a) Core and Cone 




(b) Patchy Beams 




Figure 5: Phenomenological models for pulse shape morphology produced by different 
line-of-sight cuts of the beam. (Figure designed by M. Kramer and A. von Hoensbroech) . 



Two contrasting phenomenological models to explain the observed pulse shapes are 
shown in Fig. ^. The "core and cone" model, proposed by Rankin ||205|| , depicts the beam 
as a core surrounded by a series of nested cones. Alternatively, the "patchy beam" model, 
championed by Lyne and Manchester ||165| , p.07|| , has the beam populated by a series of 
randomly-distributed emitting regions. Further work in this area, particularly in trying 
to quantify the variety of pulse shapes (number of distinct components and the relative 
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fraction that they occur) is necessary to improve our understanding of the fraction of sky 
covered by the radio pulsar emission beam. We return to this topic in the context of 
pulsar demography later on in § |3l2 . 



2.3 The pulsar distance scale 

From the sky distribution shown in Fig. |^ it is immediately apparent that pulsars are 
strongly concentrated along the Galactic plane. This indicates that pulsars populate the 
disk of our Galaxy. Unlike most other classes of astrophysical objects, quantitative esti- 
mates of the distances to each pulsar can be made from an effect known as pulse dispersion, 




Figure 6: The sky distribution of 1026 pulsars in Galactic coordinates. The plane of the 
Galaxy is the central horizontal line. The Galactic centre is the midpoint of this line. 

the delay in pulse arrival times across a finite bandwidth. Dispersion occurs because the 
group velocity of the pulsed radiation through the ionised component of the interstellar 
medium is frequency dependent: pulses emitted at higher radio frequencies travel faster 
through the interstellar medium, arriving earlier than those emitted at lower frequencies. 
The delay At in arrival times between a high frequency i^hi and a low frequency i>\o pulse 
is given [ p,72|| by 

At = 4150s X (z/i;2 _ u^^) X DM, (1) 

where the frequencies are in MHz and the dispersion measure DM (cm~^ pc) is the inte- 
grated column density of free electrons along the line of sight: 

DM = / ne dl. (2) 

Jo 

Here, d is the distance to the pulsar (pc) and rig is the free electron density (cm~'^). From 
equation |] it is obvious that a measurement of the delay across a finite bandwidth yields 
the DM. Pulsars at large distances have higher column densities and therefore larger DMs 
than those pulsars closer to Earth so that, from Equation |l|, the dispersive delay across 
the bandwidth is greater. Hence, given the DM, the distance can be estimated from a 
model of the Galactic distribution of free electrons. 

The electron density model is calibrated from the pulsars with independent distance 
estimates and measurements of scattering for lines of sight towards various Galactic and 
extragalactic sources. Independent distance estimates now exist for over 100 pulsars based 
on three basic techniques: neutral hydrogen absorption, trigonometric parallax (measured 
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either with an interferometer or through pulse time-of-arrival techniques) and from asso- 
ciations with objects of known distance {i.e. supernova remnants, globular clusters and 
the Magellanic Clouds). Based on these data, Taylor & Cordes |p40|| have developed an 
electron density model which is free from large systematic trends and can be used to 
provide distance estimates with an uncertainty of ~ 30%. However, use of this model 
to estimate distances to individual pulsars may result in uncertainties by as much as a 
factor of two. This model is currently being refined following recent pulsar discoveries and 
independent distance and scattering measurements (J. Cordes, private communication). 



2.4 Normal and millisecond pulsars 
2.4.1 Spin parameters 

The present public-domain catalogue, available on-line at Princeton |Q, contains up-to- 
date parameters for 706 pulsars. Parameters for many of the new Parkes multibeam pulsar 
surveys are also available on-line ^. Most of these are "normal" in the sense that 



their pulse periods P are of order one second and are observed to increase secularly at 
rates P of typically 10~^^ s/ s. A growing fraction of the observed sample are the so-called 
"millisecond pulsars", which have spin periods primarily in the range 1.5 and 30 ms and 
rates of slowdown < 10~^^ s/s. The first millisecond pulsar discovered, B1937-I-21 ||36| , 
with P = 1.5578 ms, remains the most rapidly rotating neutron star presently known. 
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Figure 7: The ubiquitous P — P diagram shown for a sample of radio pulsars. Those 
objects known to be members of binary systems are highlighted by a circle (for low- 
eccentricity orbits) or an ellipse (for elliptical orbits). Pulsars thought to be associated 
with supernova remnants are highlighted by the starred symbols. 



A very useful means of demonstrating the distinction between these two classes is the 
'^P P diagram" - a logarithmic scatter plot of the observed pulse period versus the period 
derivative. As shown in Fig. 0, normal pulsars occupy the majority of the upper right 
hand part of the diagram, while the millisecond pulsars reside in the lower left hand part 
of the diagram. The differences in P and P imply different ages and surface magnetic 
field strengths. By treating the pulsar as a rotating magnetic dipole, one may show that 
the surface magnetic field strength is proportional to (PP)^/^ [ [1791 ]. Lines of constant 
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magnetic field strength are drawn on Fig. together with hnes of constant characteristic 
{tc = P/C^P))- Typical inferred magnetic fields and ages are 10^^ G and 10^ yr for 
the normal pulsars and 10® G and 10^ yr for the millisecond pulsars. 



2.4.2 Binary companions 

As can be inferred from Fig. |l|, just under 4% of all known pulsars in the Galactic disk 
are members of binary systems. Timing measurements place useful constraints on 
the masses of the companions which, supplemented by observations at other wavelengths, 
tell us a great deal about their nature. The present sample of orbiting companions are 
either white dwarfs, main sequence stars, or other neutron stars. Two notable hybrid 
systems are the "planet pulsars" PSR B1257+12 and B1620-26. B1257+12 is a 6.2-ms 
pulsar accompanied by at least three Earth-mass bodies [|274 , |^, 273 1 while B1620-26, an 
11-ms pulsar in the globular cluster M4, is part of a triple system with a white dwarf and 
a high-mass planet [p49| , |3^ , p47 . 
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Figure 8: Companion mass versus orbital eccentricity for the sample of binary pulsars. 



A very important additional difference between normal and millisecond pulsars is 
the presence of an orbiting companion. Orbital companions are much more commonly 
observed around millisecond pulsars (~ 80% of the observed sample) than around the 
normal pulsars (< 1%). Fig. |^ is a scatter plot of orbital eccentricity versus mass of the 
companion. The dashed line serves merely to guide the eye in this figure. Binary systems 
lying below the line are those with low-mass companions 0.7 M0 — predominantly 
white dwarfs) and essentially circular orbits: 10~^ e < 0.1. Binary pulsars with high- 
mass companions 1 Mq — neutron stars or main sequence stars) are in eccentric orbits: 
0.15 ^ e ^ 0.9 and lie above the line. 



2.4.3 Evolutionary scenarios 

The presently favoured model to explain the formation of the various types of systems 



has been developed over the years by a number of authors |5^, ^ p25| , |22| . The model is 
sketched in Fig. ^ and is now qualitatively summarised. 

Starting with a binary star system, a neutron star is formed during the supernova 
explosion of the initially more massive star which has an inherently shorter main sequence 
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lifetime. From the virial theorem it follows that the binary system gets disrupted if more 
than half the total pre-supernova mass is ejected from the system during the explosion 



116| , 0. In addition, the fraction of surviving binaries is affected by the magnitude and 



direction of any impulsive "kick" velocity the neutron star receives at birth |116, 37 



Those binary systems that disrupt produce a high-velocity isolated neutron star and 



an OB runaway star |54]. The high binary disruption probability during the explosion 
explains, qualitatively at least, why so few normal pulsars have companions. Over the 
next 10^ yr or so after the explosion, the neutron star may be observable as a normal radio 
pulsar spinning down to a period > several seconds. After this time, the energy output of 
the star diminuishes to a point where it no longer produces significant radio emission. 

For those few binaries that remain bound, and in which the companion is sufficiently 
massive to evolve into a giant and overflow its Roche lobe, the old spun-down neutron 
star can gain a new lease of life as a pulsar by accreting matter and therefore angular 



momentum at the expense of the orbital angular momentum of the binary system |22 



The term "recycled pulsar" is often used to describe such objects. During this accretion 
phase, the X-rays produced by the liberation of gravitational energy of the infalling matter 
onto the neutron star mean that such a system is expected to be visible as an X-ray binary 
system. Two classes of X-ray binaries relevant to binary and millisecond pulsars exist, viz: 
neutron stars with high-mass or low-mass companions. For a detailed review of the X-ray 
binary population, including systems likely to contain black holes rather than neutron 
stars, the interested reader is referred to 
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Figure 9: Cartoon showing various evolutionary scenarios involving binary pulsars. 

The high-mass companions are massive enough to explode as a supernova, producing 
a second neutron star. If the binary system is lucky enough to survive this explosion, it 



ends up as a double neutron star binary. The classic example is PSR B1913-I-16 [ [L19|| , a 



8 



59-ms radio pulsar with a characteristic age of ~ 10* yr which orbits its companion every 
7.75 hr [p44| , |245|| . In this formation scenario, PSR B1913+16 is an example of the older, 
first-born, neutron star that has subsequently accreted matter from its companion. So 
far there are no clear examples of systems where the second-born neutron star is observed 
as a radio pulsar. In the case of PSR B1820— 11 | |171|| , which may be an example, the 
mass of the companion is not well determined, so either a main-sequence [ p02|| or a white 
dwarf companion [|203|| are plausible alternatives. This lack of observation of second-born 
neutron stars as radio pulsars is probably reasonable when one realises that the observable 
lifetimes of recycled pulsars are much larger than those of normal pulsars. As discussed 
in §3.4. H , double neutron star binary systems are very rare in the Galaxy — another 
indication that the majority of binary systems get disrupted when one of the components 
explodes as a supernova. Systems disrupted after the supernova of the secondary form 
a mildly-recycled isolated pulsar and a young pulsar formed during the explosion of the 
secondary. 

Although no system has so far been found in which both neutron stars are visible as 
radio pulsars, timing measurements of three systems show that the companion masses are 



1.4 Mq — as expected for neutron stars ||222|| . In addition, no optical counterparts are 



seen. Thus, we conclude that these unseen companions are neutron stars that are either 
too weak to be detected, no longer active as radio pulsars, or their emission beams do 
not intersect our line of sight. The two known young radio pulsars with main sequence 
companions massive enough to explode as a supernova probably represent the intermediate 



phase between high-mass X-ray binaries and double neutron star systems ||126| , |132 



The companions in the low-mass X-ray binaries evolve and transfer matter onto the 
neutron star on a much longer time-scale, spinning it up to periods as short as a few 
ms This model has gained strong support in recent years from the discoveries of 

quasi-periodic kHz oscillations in a number of low-mass X-ray binaries |p68| ], as well as 
Doppler-shifted 2.49-ms X-ray pulsations from the transient X-ray burster SAX J1808.4- 
3658 ||269| , [72[] . At the end of the spin-up phase, the secondary sheds its outer layers to 
become a white dwarf in orbit around a rapidly spinning millisecond pulsar. Presently 
~ 10 of these systems have compelling optical identifications of the white dwarf companion 
155| , |154|| . Perhaps the best example is the white dwarf companion to the 5.25-ms 



pulsar J1012+5307 ||187| , |150|| . This 19**^ magnitude white dwarf is bright enough to allow 
measurements of its surface gravity and orbital velocity [p59| . 



The range of white dwarf masses observed is becoming broader. Since this article 
originally appeared in 1998 the number of "intermediate-mass binary pulsars" has 
grown significantly These systems are distinct to the "classical" millisecond pulsar- 
white dwarf binaries like PSR J1012+5307 in several ways: (1) the spin period of the 
radio pulsar is generally longer (9-200 ms); (2) the mass of the white dwarf is larger 
(typically close to 1 Mq); (3) the orbital eccentricity, while still essentially circular, is often 
significantly larger (~ 10^^). It is not presently clear whether these systems originated 
from either low- or high-mass X-ray binaries. It was suggested by van den Heuvel 



that they have more in common with high-mass systems, the difference being that the 
secondary star was not sufficiently massive to explode as a supernova. Instead it formed 
a white dwarf. Detailed studies of this sub-population of binary pulsars are required for 
further understanding in this area. 

Another relatively poorly understood area is the existence of solitary millisecond pul- 
sars in the Galactic disk (which comprise just under 20% of all Galactic millisecond 
pulsars) . Although it has been proposed that the millisecond pulsars have got rid of their 



companion by ablation, as appears to be happening in the PSR B1957-I-20 system |101 
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it is not clear whether the time-scales for this process are feasible. There is some obser- 
vational evidence that suggests that solitary millisecond pulsars are less luminous than 
binary millisecond pulsars |139|| . If confirmed by future discoveries, this would need 
to be explained by any viable evolutionary model. 



2.4.4 Space velocities 

Pulsars have long been known to have space velocities at least an order of magnitude larger 
than those of their main sequence progenitors, which have typical values between 10 and 
50 km s~^. The first direct evidence for large velocities came from optical observations 



of the Crab pulsar in 1968 ||253|| , showing that the neutron star has a velocity in excess 
of 100 km s~^. Proper motions for about 100 pulsars have subsequently been measured 
largely by radio interferometric techniques [|15^ , 40, 96, |110|] . These data imply a broad 



velocity spectrum ranging from to over 1000 km s [p.64|| 




Figure 10: A simulation following the motion of 100 pulsars in a model gravitational 
potential of our Galaxy for 200 Myr. The view is edge-on i. e. the horizontal axis represents 
the Galactic plane (30 kpc across) while the vertical axis represents ±10 kpc from the 
plane. This snapshot shows the initial configuration of young neutron stars. Click here 
to see the movie in action. 

Such large velocities are perhaps not surprising, given the violent conditions under 
which neutron stars are formed. Shklovskii [p24|| demonstrated that, if the explosion is 
only slightly asymmetric, an impulsive "kick" velocity of up to 1000 km s~^ is imparted 
to the neutron star. In addition, if the neutron star progenitor was a member of a binary 
system prior to the explosion, the pre-supernova orbital velocity will also contribute to 
the resulting speed of the newly-formed pulsar. High-velocity pulsars born close to the 
Galactic plane quickly migrate to higher Galactic latitudes. This migration is seen in 
Fig. a dynamical simulation of the orbits of 100 neutron stars in a model of the 
Galactic gravitational potential. Given such a broad velocity spectrum, as much as half 
of all pulsars will eventually escape the gravitational potential of the Galaxy and end up 
in intergalactic space ||164| , [f7[] . 

Based on the proper motion data, recent studies have demonstrated that the mean 
birth velocity of normal pulsars is ~ 450 km s~^ ( ||164| , |149| , [77| , p.02|| ; see, however, also 
111| , |109|] ). This is significantly larger than the velocities of millisecond and binary pulsars. 
Recent studies suggest that their mean birth velocity is likely to be in the range ~ 80 — 140 



km s [145, 76, 168]. The main reason for this difference surely lies in the fact that about 
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80% of the millisecond pulsars are members of binary systems (§ |2.4| ) which could not have 
survived had the neutron star received a substantial kick velocity. 



2.5 Searching for pulsars 

Pulsar searching is, conceptually at least, a rather simple process — the detection of 
a dispersed, periodic signal hidden in a noisy time series collected using a large radio 
telescope. The search These data are then analysed for periodic signals. We give here a 
brief description of the basic search techniques. More detailed discussions can be found 
elsewhere jlSTl, [l88i ra . 



A schematic pulsar search is shown in Fig. 11. The finite bandwidth is split up into 
a number of channels, typically using a filterbank or a correlator (see e.g. |3^), either 
of which usually provides a much finer frequency channelisation than the eight channels 
shown for illustrative purposes in Fig. |ll[ The channels are then de-dispersed (see § p.l.2|) 
to form a single noisy time series. An efficient way to find a periodic signal in these data is 
to take the Fast Fourier Transform (FFT) and plot the resulting amplitude spectrum. For 
a narrow pulse the spectrum will show a family of harmonics. To detect weaker signals 
still, a harmonic summing technique is usually implemented ||157|| . The best candidates 
are saved and the whole process is repeated for another trial DM. 
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Figure 11: Schematic summarising the essential steps in a "standard" pulsar search. 



After the data have been processed for a suitable range of DM, a list of pulsar candi- 
dates is compiled and the raw time series data are folded modulo each candidate period. 
In practise the analysis is often hampered by the presence of periodic interference sources 
which can often look very "pulsar-like" . Although interference excision schemes (usually 
based on coincidence analyses of data taken from different points on the sky) work fairly 
well, interference is an ever-increasing problem in radio astronomy and considerable efforts 
are required to carry out sensitive searches. 



2.6 Where to look for binary and millisecond pulsars 

Following the above discussions on demography and evolution, it is instructive to briefly 
summarize the rationale behind the major searches being carried out at the present time. 



11 



2.6.1 All-sky searches 

The oldest radio pulsars form a virialised population of stars oscillating in the Galactic 
gravitational potential. The scale height for such a population is at least 500 pc, about 
10 times that of the massive stars which populate the Galactic plane. Since the typical 
ages of millisecond pulsars are several Gyr or more, we expect, from our vantage point in 
the Galaxy, to be in the middle of an essentially isotropic population of nearby sources. 
All-sky searches for millisecond pulsars at high Galactic latitudes have been very effective 
in probing this population. Much of the initial interest and excitement in this area 
was started at Arecibo when Wolszczan discovered two classic recycled pulsars at high 
latitudes: the neutron star binary B 1534+ 12 [p71|] and the planets pulsar B 1257+ 12 ||274 



Surveys carried out at Arecibo, Parkes, Jodrell Bank and Green Bank by others in the 
1990s found many other millisecond and recycled pulsars in this way. Camilo has written 
several excellent reviews of these surveys |^ |6^, 0. See also Tables 0, |^ and || in the 
appendix. 



2.6.2 Searches close to the plane of our Galaxy 

Young pulsars are most likely to be found near to their place of birth, close to the Galactic 
plane. This is the target region of one of the Parkes multibeam surveys and has already 
resulted in the discovery of around 600 new pulsars |175|| , almost half the number 



currently known! Such a large haul inevitably results in a number of interesting individual 
objects such as: PSR J1141— 6545, a young pulsar in a relativistic 4-hr orbit around a 
white dwarf | |133|| ; PSR J1740— 3052, a young pulsar orbiting an ~ 11 Mq star (probably 
a giant ||175|| ); several intermediate-mass binary pulsars ||68| and a likely double neutron 
star system PSR J1811-1736 Wm. 



2.6.3 Searches at intermediate Galactic latitudes 

In order to probe more deeply into the population of millisecond and recycled pulsars 
than possible at high Galactic latitudes, Edwards et al. have recently completed a survey 
of intermediate latitudes with the Parkes multibeam system ^2, 93]. The results of this 
survey are extremely exciting — 58 new pulsars including 8 relatively distant recycled 
objects. Two of the new recycled pulsars from this survey recently announced by Edwards 
& Bailes 
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are mildly relativistic neutron star-white dwarf binaries. An analysis of the 
full results from this survey should significantly improve our knowledge on the Galaxy- 
wide population and birth-rate of millisecond pulsars. 



2.6.4 Targeted searches of globular clusters 

Globular clusters have long been known to be breeding grounds for millisecond and binary 
pulsars. The main reason for this is the high stellar density in globular clusters relative to 
most of the rest of the Galaxy. As a result, low-mass X-ray binaries are almost 10 times 
more abundant in clusters than in the Galactic disk. In addition, exchange interactions 
between binary and multiple systems in the cluster can result in the formation of exotic 
binary systems. Since a single globular cluster usually fits well within a single telescope 
beam, deep targeted searches can be made. Once the DM of a pulsar is known in a 
globular cluster, the DM parameter space for subsequent searches is essentially fixed. 
This allows computation power to be invested in so-called acceleration searches for short- 



period binary systems (see § |3.1.3| ). To date, searches have revealed 47 pulsars in globular 
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clusters (see Table § in the appendix for a list and the review by Kulkarni & Anderson 
|140|| ). Highlights include the double neutron star binary in M15 |p04|| and a low- mass 
binary system with a 95-min orbital period in 47 Tucanae 
pulsars currently known in this cluster alone |^ 



66 , one of 20 millisecond 



On-going surveys of clusters continue 



to yield new discoveries [206, 81 



2.7 Going further 



Two excellent graduate-level monographs are available: the classic 1970s text Pulsars 
by Manchester & Taylor [|179|| and the more up-to-date Pulsar Astronomy by Lyne & 
Smith | 172 |. Those wishing to approach the subject from a more theoretical viewpoint 
are advised to read Michel's The Theory of Neutron Star Magnetospheres ||180|] and The 
Physics of the Pulsar Magnetosphere by Beskin, Gurevich & Istomin . Our summary 
of evolutionary aspects serves merely as a primer to the vast body of literature available. 
The reader is referred to the excellent review by Bhattacharya & van den Heuvel 
for further insights. For an excellent overview of pulsar distance measurements and their 
implications, see the review by Weisberg [p65|] . 

Pulsar resources available on the Internet are continually becoming more extensive 
and useful. Good starting points for pulsar- surfers are the pages maintained at Arecibo 
@], Berkeley [§, Bonn g, Jodrell Bank 0, Princeton 



Swinburne [10] and Sydney 
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3 The Galactic Pulsar Population 



Soon after the discovery of pulsars, it was realised that the observed sample is heavily 
biased towards the brighter objects that are the easiest to detect. What we observe 
therefore most likely represents only the tip of the iceberg of a much larger underlying 
population ||106|| . The extent to which the sample is incomplete is well demonstrated by 
the projection of pulsars onto the Galactic plane and their cumulative number distribution 
as a function of distance shown in Fig. O. Although the clustering of sources around the 



Sun seen in the left panel of Fig. |I2| would be consistent with Ptolemy's geocentric picture 
of the heavens, it is clearly at variance with what we now know about the Galaxy, where 
the massive stars show a radial distribution about the Galactic center. 
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Figure 12: Left: The sample of radio pulsars from the Princeton catalog Q projected 
onto the Galactic plane. The Galactic center is at (0,0) and the Sun is at (-8.5,0). Right: 
Cumulative number of observed pulsars (solid line) as a function of projected distance, d. 
The dashed line shows the expected distribution for a model population (see text). 



The extent to which the pulsar sample is incomplete is shown in the right panel of 
Fig. O where the cumulative number of pulsars is plotted as a function of the projected 
distance from the Sun. The observed distribution is compared to the expected distribution 
for a simple model population in which there are errors in the distance scale, but no 
selection effects. We see that the observed sample becomes strongly deficient in terms of 
the number of sources for distances beyond a few kpc. We now discuss the main selection 
effects that hamper the detection of pulsars in some detail. 



3.1 Selection effects in pulsar searches 
3.1.1 The inverse square law and survey thresholds 

The most prominent selection effect at play in the observed pulsar sample is the inverse 
square law, i.e. for a given intrinsic luminosity^, the observed flux density falls off as the 
inverse square of the distance. This results in the observed sample being dominated by 
nearby and/or bright objects. Beyond distances of a few kpc from the Sun, the apparent 

^Pulsar astronomers usually define the luminosity L = S(P, where S is the mean flu x density at 400 
MHz (a standard observing frequency) and d is the distance derived from the DM (§2.2). 
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flux density falls below the flux thresholds Smin of most surveys. Following [B£], we write: 

-1 - A . -Ml . -1/2 ^ ^ ^1/2 



5n 



SNR, 



Tree ~l~ -^f 



sky 



K 



G 



KJy-\ 



VMHz/ 



P-W 



mJy, 



(3) 

where SNRmin is the threshold signal-to-noise ratio, Tree and T^^y are the receiver and sky 
noise temperatures, G is the gain of the antenna, Ai/ is the observing bandwidth, tint is 
the integration time, W is the detected pulse width and P is the pulse period. 



3.1.2 Pulse dispersion and scattering 

It follows from Equation |] that the sensitivity decreases as W/ {P — W) increases. Also 
note that if W ^ P, the pulsed signal is smeared into the background emission and is 
no longer detectable, regardless of how luminous the source may be. The detected pulse 
width W will be broader than the intrinsic value largely as a result of pulse dispersion and 
scattering by free electrons in the interstellar medium. As discussed above, the dispersive 
smearing scales as Az//z/^, where u is the observing frequency. This can largely be removed 
by dividing the pass-band into a number of channels and applying successively longer time 
delays to higher frequency channels before summing over all channels to produce a sharp 
profile. This process is known as incoherent dedispersion. 

The smearing across the individual frequency channels, however, still remains and be- 
comes significant at high dispersions when searching for short-period pulsars. Multi-path 
scattering results in a one-sided broadening due to the delay in arrival times which scales 
roughly as which can not be removed by instrumental means. A simple scattering 
model is shown in Fig. |l^ in which the scattering electrons are assumed to lie in a thin 



screen between the pulsar and the observer ||217 | 
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Figure 13: Pulse scattering caused by irregularities in the interstellar medium. The 
difference in path lengths and therefore in arrival times of the scattered rays results in a 
^'scattering tail" in the observed pulse profile which lowers its signal-to-noise ratio. 

Dispersion and scattering are most severe for distant pulsars in the inner Galaxy where 
the number of free electrons along the line of sight becomes large. The strong frequency 
dependence of both effects means that they are considerably less of a problem for surveys 



at observing frequencies > 1400 MHz |7^, |125|| compared to the usual 400 MHz search 
frequency. An added bonus for such observations is the reduction in Tgky, since the spectral 
index of the non-thermal Galactic emission is about -2.8 ||143|| . Pulsars themselves have 
steep radio spectra. Typical spectral indices are -1.6 



so that flux densities are 
roughly an order of magnitude lower at 1400 MHz compared to 400 MHz. Fortunately, 
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this can usually be compensated for by the use of larger receiver bandwidths at higher 
radio frequencies. For example, the 1370-MHz system at Parkes has a bandwidth of 288 
MHz [ |163| ] compared to the 430-MHz system, where 32 MHz is available [ |176|| . 



3.1.3 Orbital acceleration 

Standard pulsar searches use Fourier techniques to search for a-priori unknown periodic 
signals and usually assume that the apparent pulse period remains constant throughout 
the observation. For searches with integration times much greater than a few minutes 
this assumption is only valid for solitary pulsars, or those in binary systems where the 
orbital periods are longer than about a day. For shorter-period binary systems, as noted 
by Johnston & Kulkarni 



the Doppler-shifting of the period results in a spreading 
of the signal power over a number of frequency bins in the Fourier domain, leading to a 
reduction in signal-to-noise ratio. An observer will perceive the frequency of a pulsar to 
shift by an amount aT/{Pc), where a is the (assumed constant) line-of-sight acceleration 
during the observation of length T, P is the (constant) pulse period in its rest frame and 
c is the speed of light. Given that the width of a frequency bin is 1/T, we see that the 
signal will drift into more than one spectral bin if aT'^/{Pc) > 1. Survey sensitivities to 
rapidly-spinning pulsars in tight orbits are therefore significantly compromised when the 
integration times are large. 
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Figure 14: Left: a 22.5-min Arecibo observation of the binary pulsar B1913-hl6. The 
assumption that the pulsar has a constant period during this time is clearly inappropriate 
given the drifting in phase of the pulse during the integration (grey scale plot). Right: the 
same observation after applying an acceleration search. This shows the effective recovery 
of the pulse shape and a significant improvement in the signal-to-noise ratio. 



As an example of this effect, as seen in the time domain. Fig. |IJ shows a 22.5-min search 
mode observation of Hulse & Taylor's famous binary pulsar B1913+16 ||119| , P44| , |245 . 
Although this observation covers only about 5% of the orbit (7.75 hr), the effects of the 
Doppler smearing on the pulse signal are very apparent. While the standard search code 
(seeking constant periodicity) nominally detects the pulsar with a signal-to-noise ratio of 
9.5 for this observation, it is clear that the Doppler shifting of the pulse period seen in 
the individual sub-integrations results in a significant reduction in signal-to-noise. 

It is clearly desirable to employ a technique to recover the loss in sensitivity due to 



Doppler smearing. One such technique, the so-called "acceleration search" [|181| 
the pulsar has a constant acceleration during the integration. Each time series can then 
be re-sampled to refer it to the frame of an inertial observer using the Doppler formula 
to relate a time interval, r, in the pulsar frame to that in the observed frame at time t, 



16 



as r(t) oc (1 + at/c). Searching over a range of accelerations is desirable to find the time 
series for which the trial acceleration most closely matches the true value. In the ideal 
case, a time series is produced with a signal of constant period for which full sensitivity 
is recovered (see right panel of Fig. |I^). Anderson et al. ||2^ used this technique to find 



PSR B2127+11C, a double neutron star binary in M15 which has parameters similar to 
B1913+16. Camilo et al. |66] have recently applied the same technique to 47 Tucanae to 



discover 9 binary pulsars, including one in a 96-min orbit around a low-mass (0.15 Mq) 
companion. This is currently the shortest binary period for any known radio pulsar. 

For the shortest orbital periods, the assumption of a constant acceleration during 
the observation clearly breaks down. Ransom et al. [|206|| have developed a particularly 



efficient algorithm for finding binaries whose orbits are so short that many orbits can take 
place during an integration. This phase modulation technique exploits the fact that the 
pulses are modulated by the orbit to create a family of periodic sidebands around the 
nominal spin period of the pulsar. This technique has already been used to discover a 
1.7-hr binary pulsar in NGC 6544 ||206|| . The existence of these short-period radio pulsar 



binaries, as well as the 11-min X-ray binary X1820— 303 in NGC 6624 ||233|| , implies 



that there must be many more short-period binaries containing radio or X-ray pulsars in 
globular clusters that are waiting to be discovered by more sensitive searches. 

3.2 Correcting the observed pulsar sample 
3.2.1 Scale Factor determination 

Now that we have a fiavour for the variety and severity of the selection effects that plague 
the observed sample of pulsars, how do we decouple these effects to form a less biased 
picture of the true population of objects? A very useful technique, first employed by 
Phinney & Blandford and Vivekanand & Narayan ||200| , |261|| , is to define a scaling factor 



^ as the ratio of the total Galactic volume weighted by pulsar density to the volume in 
which a pulsar could be detected by the surveys: 

' ' JJpL T.{R,z)RdRdz ■ ^ ' 

In this expression, E(i?, z) is the assumed pulsar distribution in terms of galactocentric 
radius R and height above the Galactic plane z. Note that is primarily a function of 
period P and luminosity L such that short period/low-luminosity pulsars have smaller 
detectable volumes and therefore higher ^ values than their long period/high-luminosity 
counterparts. This approach is similar to the classic V^/l^ax technique first used to correct 
observationally-biased samples of quasars [|218 |. 



This technique can be used to estimate the total number of active pulsars in the 
Galaxy. In practice, this is achieved by calculating ^ for each pulsar separately using a 
Monte Carlo simulation to model the volume of the Galaxy probed by the major surveys 



183|. For a sample of Ai'obs observed pulsars above a minimum luminosity Lmin, the total 



number of pulsars in the Galaxy with luminosities above this value is simply 

E J, (5) 

i=l 



where / is the model-dependent "beaming fraction" discussed below in §3.2.3. Monte 



Carlo simulations of the pulsar population incorporating the aforementioned selection 



effects have shown this method to be reliable, as long as iVobs is reasonably large [148|. 
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3.2.2 The small-number bias 



For small samples of observationally-selected objects, the detected sources are likely to be 
those with larger-than-average luminosities. The sum of the scale factors (Equation 
therefore, will tend to underestimate the true size of the population. This "small-number 
bias" was first pointed out by Kalogera et al. ||129| , |130|| for the sample of double neutron 
star binaries where we know of only three clear-cut examples (§ pj.4.1| ). Only when the 
number of sources in the sample gets past 10 or so does the sum of the scale factors 
become a good indicator of the true population size. 



Model Galactic Population 




Figure 15: Small-number bias of the scale factor estimates derived from a synthetic pop- 
ulation of sources where the true number of sources is known. Left: an edge-on view of 
a model Galactic source population. Right: the thick line shows Nq, the true number of 
objects in the model Galaxy, plotted against N^hs, the number detected by a fiux-limited 
survey. The thin solid line shows N^st, the median sum of the scale factors, as a func- 
tion of No\^s from a large number of Monte- Carlo trials. Dashed lines show 25 and 75% 
percentiles of the N^st distribution. 



3.2.3 The beaming correction 

The "beaming fraction" / in Equation ^ is simply the fraction of An steradians swept out 
by the radio beam during one rotation. Thus / gives the probability that the beam cuts 
the line-of-sight of an arbitrarily positioned observer. A naive estimate of / is 20%; this 
assumes a beam width of ~ 10° and a randomly distributed inclination angle between 



the spin and magnetic axes [p43|| . Observational evidence suggests that shorter period 
pulsars have wider beams and therefore larger beaming fractions than their long-period 
counterparts ||184| , |165| , |^ , [^36|| . It must be said, however, that a consensus on the 



beaming fraction-period relation has yet to be reached. This is shown in Fig. ^ where 
we compare the period dependence of / as given by a number of models. Adopting the 
Lyne & Manchester model, pulsars with periods ~ 100 ms beam to about 30% of the sky 
compared to the Narayan & Vivekanand model in which pulsars with periods below 100 
ms beam to the entire sky. 

When most of these beaming models were originally proposed, the sample of millisec- 
ond pulsars was ^ 5 and hence their predictions about the beaming fractions of short- 
period pulsars relied largely on extrapolations from the normal pulsars. A recent analysis 
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Figure 16: Beaming fraction plotted against pulse period for four different beaming mod- 
els: Tauris & Manchester 1998 (TM88; pWj] ), Lyne & Manchester 1988 (LM88; 
Biggs 1990 (JDB90; 1^) and Narayan & Vivekanand 1983 (NV83 Jl8i|/ ). 



of a large sample of millisecond pulsar profiles by Kramer et al. [|139|| suggests that the 
beaming fraction of millisecond pulsars lies between 50 and 100%. 



3.3 The population of normal and millisecond pulsars 

3.3.1 Luminosity distributions and local number estimates 

The most recent use of the scale factor approach to derive the characteristics of the true 
normal and millisecond pulsar populations is based on the sample of pulsars within 1.5 kpc 
of the Sun 
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The rationale for this cut-off is that, within this region, the selection 
effects are well understood and easier to quantify by comparison with the rest of the 
Galaxy. These calculations should give reliable estimates for the local pulsar population. 





10 100 
Luminosity (mJy kpc^) 



1000 



1 10 

Luminosity (mJy kpc^) 



100 



Figure 17: Left: The corrected luminosity distribution (solid histogram with error bars) 
for normal pulsars. The corrected distribution before the beaming model has been applied 
is shown by the dot-dashed line. Right: The corresponding distribution for millisecond 
pulsars. In both cases, the observed distribution is shown by the dashed line and the 
thick solid line is a power law with a slope of -1. The difference between the observed and 
corrected distributions highlights the severe under- sampling of low- luminosity pulsars. 
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The luminosity distributions obtained from this analysis are shown in Fig. Tj. For 
the normal pulsars, integrating the corrected distribution above 1 mJy kpc^ and dividing 



by vr X (1.5)^ kpc^ yields a local surface density, assuming Biggs' beaming model of 
156 ± 31 pulsars kpc~^. The same analysis for the millisecond pulsars, assuming a mean 



beaming fraction of 75% ||139|| , leads to a local surface density of 38 ± 16 pulsars kpc ^ 
for luminosities above 1 mJy kpc^. 

3.3.2 Galactic population and birth-rates 

Integrating the local surface densities of pulsars over the whole Galaxy requires a knowl- 



edge of the presently rather uncertain Galactocentric radial distribution ||169| , |123|| . One 
approach is to assume that pulsars have a radial distribution similar to that of other 
stellar populations and scale the local number density with this distribution to estimate 
the total Galactic population. The corresponding local-to-Galactic scaling is 1000±250 
kpc^ [|207|| . With this approach we estimate there to be ~ 160,000 active normal pulsars 



and ~ 40, 000 millisecond pulsars in the Galaxy. Based on these estimates, we are in 
a position to deduce the corresponding rate of formation or birth-rate. From the P P 
diagram in Fig. |^, we infer a typical lifetime for normal pulsars of ~ 10^ yr, correspond- 
ing to a Galactic birth rate of ~ 1 per 60 yr — consistent with the rate of supernovae 
||255|| . As noted in § p^.4| , the millisecond pulsars are much older, with ages close to that 
of the Universe, Tu (we assume here r„ = 10 Gyr ||121|| ). Taking the maximum age of the 



millisecond pulsars to be r„, we infer a mean birth rate of at least 4 x 10 yr . This is 



consistent, within the uncertainties, with the birth-rate of low-mass X-ray binaries 152 



3.3.3 Implications for gravitational wave detectors 

The estimates of the local surface density of active pulsars allow us to deduce the likely 
distance to the nearest neutron star to Earth. For the combined millisecond and nor- 
mal pulsar populations, with a surface density of 193 ± 35 pulsars kpc~^, the nearest 
neutron star is thus likely to be ^ 40 pc. This number is of interest to those build- 
ing gravitational wave detectors, since it determines the likely amplitude of gravitational 
waves emitted from nearby rotating neutron stars p20| |. According to Thorne ||246|| , cur- 



rently planned detectors will be able to detect neutron stars with ellipticities greater than 
7.5 X 10~^^(P(i)^, where P is the rotation period in ms and d is the distance in kpc. 



The recent probable detection of free precession in the radio pulsar B 1828— 11 ||228|| does 
indicate that ellipticities exist in neutron stars so that nearby objects may be continuous 
sources of gravitational radiation. 

Thus, in order to detect small ellipticities, nearby sources with short spin periods are 
required. One of the best candidates is the nearby 5.75-ms pulsar J0437-4715 [T^]. At a 



distance of 178 it 26 pc |pi5|| this is currently the closest known millisecond pulsar to the 



Earth. The closest known neutron star is RX J185635— 3754 discovered in the ROSAT 



all-sky survey p63|| . Multi-epoch HST observations show that this isolated neutron star is 
located at a distance of 61 ±9 pc [p62|| . In keeping with other radio-quiet isolated neutron 



stars, the period of this pulsar is likely to be several seconds | 186 |. 



3.4 The population of relativistic binaries 

Of particular interest to the astronomical community at large are the numbers of relativis- 
tic binary systems in the Galaxy. Systems involving neutron stars are: double neutron 
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star binaries, neutron star-white dwarf binaries and neutron star-black hole binaries. In- 
terest in these systems is two-fold: (1) to test the predictions of general relativity against 
alternative theories of strong-field gravity using the radio pulsar as a highly stable clock 
moving in the strong gravitational field; (2) to detect strong gravitational wave emission 
from coalescing binaries with upcoming gravitational-wave observatories like GEO600 and 



LIGO, VIRGO and TAMApl6|, pT9 



Although no radio pulsar in orbit around a black hole companion has so far been 
observed, we now know of several double neutron star and neutron-star white dwarf 
binaries which will merge due to gravitational wave emission within a reasonable time- 
scale. The merging time Tg of a binary system containing two compact objects due to the 
emission of gravitational radiation can be calculated from the following formula which 
requires only the component masses and current orbital period and eccentricity e: 

Here mi^2 are the masses of the two stars and /i = mim2/(mi + 1112). This formula is a 
good analytic approximation (within a few percent) to the numerical solution of the exact 
equations for Tg in the original papers by Peters & Mathews ||196| , |197|| . In the following 



subsections we review current knowledge on the population sizes and merging rates of 
such binaries where one component is visible as a radio pulsar. 

3.4.1 Double neutron star binaries 



As noted in §2^, double neutron star (DNS) binaries are expected to be rare since the 
binary system has to survive two supernova explosions. This expectation is certainly borne 
out by radio pulsar searches which have revealed only three certain DNS binaries so far: 



PSRs B1534+12 127T|, B1913+16 [ig and B2127+11C [gg. Although we cannot see 



the companion neutron star in any of these systems, we are "certain" of the identification 
from the precise measurements of the component masses via relativistic effects measured 
in pulsar timing observations (see § |4.1| ). The spin and orbital parameters of these pulsars 
are listed in Table |l|. Also listed in Table are three further DNS candidates with eccentric 



J1518+4904 B1534+12 J1811-1736 B1820-11 B1913+16 B2127+11C 



P (ms) 


40.9 


37.9 


104.2 


279.8 


59.0 


30.5 


n (d) 


8.6 


0.4 


18.8 


357.8 


0.3 


0.3 




0.25 


0.27 


0.83 


0.79 


0.62 


0.68 


Tc (10« y) 


200 


2.5 


970 


0.04 


1.1 


0.97 


r, (10« y) 


» T, 


27 


> T„ 


> T„ 


3.0 


2.2 


Masses measured? 


No 


Yes 


No 


No 


Yes 


Yes 



Table 1: Known DNS binaries and candidates. Listed are the pulse period P, the orbital 
period Pb, the orbital eccentricity e, the pulsar characteristic age t^, the expected binary 
coalescence time-scale Tg due to gravitational wave emission calculated from Equation 
Cases for which Tg is a factor of 100 or more greater than the age of the Universe are 
listed as ^ Tu- To distinguish between definite and candidate DNS systems, we also list 
whether the masses of both components have been determined. 

orbits and large mass functions but for which there is presently not sufficient component 
mass information to confirm their nature. 
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Despite the uncertainties in identifying DNS binaries, for the purposes of determining 
the Galactic merger rate, the systems for which Tg is less than Tu (i.e. PSRs B1534+12, 
B1913+16 and B2127+11C) are primarily of interest. Of these PSR B2127+11C is in the 
process of being ejected from the globular cluster M15 ||204| P01|| and is thought to make 
only a negligible contribution to the merger rate [|199|| . The general approach with the 



remaining two systems is to derive scale factors for each object (as outlined in § p.2.1 ) and 
then divide these by a reasonable estimate for the lifetime. In what follows we summarize 
the main studies of this kind. The most comprehensive investigation of the DNS binary 
population to date is the recent study by Kalogera et al. (hereafter KNST; [|13CI|| ). 

As discussed in § |3.2.1| , scale factors are dependent on the assumed pulsar distribution. 
The key parameter here is the scale height of the population with respect to the Galactic 
plane which itself is a function of the velocity distribution of the population. KNST 
examined this dependence in detail and found scale heights in the range 0.8-1.7 kpc. 



Based on this range, KNST revised earlier scale factor estimates [|79[ to 145-200 for 
B1534+12 and 45-60 for PSR B1913+16. As mentioned in § pX^ scale factors calculated 
from a small sample of objects are subject to a significant bias. KNST find the bias in 
their sample to be anywhere between 2 and 200. This boosts the scale factors to the range 
190-40000 for B1534+12 and 90-12000 for B1913+16. 

The above scale factors also require a beaming correction. As noted in § |3.2.3| , current 
radio pulsar beaming models vary considerably. Fortunately, for the two pulsars under 
consideration, detailed studies of the beam sizes H, |13|, m lead KNST to conclude 



that both pulsars beam to only about a sixth of the entire sky. The beaming-corrected 
numbers suggest a total of between 1680 and 312,000 active DNS binaries in our Galaxy. 
Many of these systems will be extremely faint objects. These estimates are dominated by 
the small-number bias factor. KNST's study highlights the importance of this effect. 
Some debate exists about what is the most reasonable estimate of the lifetime. Phin- 



ney ||199|| defines this as the sum of the pulsar's spin-down age plus Tg defined above. 
A few years later, van den Heuvel and myself argued | 257 | that a more likely estimate 
can be obtained by appealing to steady-state arguments where we expect sources to be 
created at the same rate at which they are merging. The mean lifetime was then found to 
be about three times the current spin-down age. This argument does, however, depend 
on the luminosity evolution of radio pulsars which is currently only poorly understood. 
Arzoumanian, Cordes & Wasserman |^| used kinematic data to constrain the most likely 
ages of the DNS binaries. They note that the remaining detectable lifetime should also 
take account of the reduced detectability at later epochs due to acceleration smearing 
as the DNS binary becomes more compact due to gravitational wave emission. KNST 
concluded that the lifetimes are dominated by the latter time-scale which, following Ar- 
zoumanian et al., they took to be the time for the orbital period to halve. The resulting 
lifetimes are 2.5 x 10^ yr for B1534+12 and 2.5 x 10^ yr for B1913+16. 

Taking these number and lifetime estimates, KNST find the Galactic merger rate of 
DNS binaries to range between 3 x 10~^ and 4 x 10~^ yr~^. Extrapolating this number out 
to include DNS binaries detectable by LIGO in other galaxies a la Phinney [|199|| KNST 
find the expected event rate to be < 0.25 yr~^ for LIGO-I and 2-1300 yr~^ for LIGO-II. 
Thus, despite the uncertainties, it seems that the prospects for detecting gravitational- 
wave emission from DNS inspirals in the near future are most promising. 
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3.4.2 White dwarf— neutron star binaries 



The population of white dwarf-neutron star (WDNS) sources containing a young radio 
pulsar has only recently been confirmed by observers following the identification of a white 
dwarf companion to the binary pulsar B2303+46 by van Kerkwijk et al. [p58|| . Previously, 
this eccentric binary pulsar was thought to be an example of a DNS binary in which the 
visible pulsar is the second-born neutron star | f^35| , |159|| . The optical identification rules 
this out and now strongly suggests a scenario in which the white dwarf was formed first. 
In this case, material was transfered onto the secondary during the giant phase of the 
primary so that the secondary became massive enough to form a neutron star ||258| , |237 . 



PSR B2303+46 has a long orbital period and does not contribute significantly to the 
overall merger rate of WDNS binaries. The new discovery of PSR J1141-6545 ||133|| , 
which will merge due to gravitational- wave emission within 1.3 Gyr, is suggestive of a 
large population of similar binaries. This is particularly compelling when one considers 
that the radio lifetime of the visible pulsar is only a fraction of total lifetime of the binary 
before coalescence due to gravitational-wave emission. Edwards & Bailes [^] estimate 
there to be 850 WDNS binaries within 3 kpc of the Sun which will merge within r^. 

Population syntheses by Tauris & Sennels [|237|| suggest that the formation rate of 



WDNS binaries is between 10-20 times that of DNS binaries. Based on the merging rate 
estimates for DNS binaries discussed in the previous section, this translates to a merging 
rate of WDNS binaries of between 3 x 10~^ yr~^ and 8 x 10"'^ yr~^- In summary, although 
statistics are necessarily poor at this stage, coalescing WDNS binaries look to be very 
promising sources for gravitational wave detectors. 

3.5 Going further 

Studies of pulsar population statistics represent a large proportion of the pulsar literature. 
During this section we have tried to cite many of the key papers in this field. Good starting 



points for further reading can be found in other review articles ||180|, pOj , |172|| . Our coverage 
of compact object coalescence rates has concentrated on empirical methods and we have 
hopefully convinced the reader that these are fair and straightforward. 

An alternative approach is to undertake a full-blown Monte Carlo simulation of the 
most likely evolutionary scenarios described in § |2.4.3| . In this "scenario-machine" ap- 
proach, a population of primordial binaries is synthesized with a number of underlying 
distribution functions: primary mass, binary mass ratio, orbital period distribution etc. 
The evolution of both stars is then followed to give a predicted sample of binary systems of 
all the various types. Since the full range of binary parameters is known, the merger rates 
of each type of binary are then automatically predicted by this model without the need to 
debate what the likely coalescence times will be. Selection effects are not normally taken 
into account in this approach. The final census is usually normalized to the star formation 
rate. Numerous examples of the scenario-machine approach (most often to populations of 



binaries where one or both members are NSs) can be found in the literature |211| , |254] 



These include the widely-cited code, developed by Lipunov and collaborators to perform 
population syntheses of binary stars |jll|. Although extremely instructive, the uncertain 
assumptions about initial conditions, the physics of mass transfer and the kicks applied 
to the compact object at birth result in a wide range of predicted event rates which are 
currently broader than the empirical methods [ |129|| . Ultimately, the detection statistics 
from the gravitational wave detectors could provide far tighter constraints on the DNS 
merging rate than the pulsar surveys from which these predictions are made. 
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An excellent overview of gravitational-wave astronomy and the detection of gravita- 
tional waves from inspiraling binaries is presented by Thorne in his presentation at the 



centennial meeting of the American Physical Society which is available on-line |T2 
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4 Pulsar Timing 



It became clear soon after their discovery that pulsars are excellent celestial clocks. 

TTSj 



the original discovery paper 



the period of the first pulsar to be discovered, 



In 
PSR 



B1919+21, was found to be stable to one part in 10^ over a time-scale of a few months. 
Following the discovery of the millisecond pulsar B1937-I-21 in 1982 it was demon- 
strated that its period could be measured to one part in 10^^ or better |8^. This unrivaled 
stability leads to a host of applications including time keeping, probes of relativistic grav- 
ity and natural gravitational wave detectors. 



4.1 Observing basics 

As each new pulsar is discovered, the standard practice is to add it to a list of pulsars 
which are regularly observed at least once or twice per month by large radio telescopes 
throughout the world. The schematic diagram shown in Fig. |18| summarises the essential 
steps involved in such a "time-of-arrival" (TOA) measurement. Incoming pulses emitted 
by the rotating neutron star traverse the interstellar medium before being received by 
the radio telescope. After amplification by high sensitivity receivers, the pulses are de- 
dispersed and added to form a mean pulse profile. 



Rotation axis 

, Mean Pulse Profile 
- — ^ ~ — ^ Reference clock , ^ 




Receiver ^ ' 

Figure 18: Schematic showing the main stages involved in pulsar timing observations. 



During the observation, the data regularly receive a time stamp, usually based on 
a caesium time standard or hydrogen maser at the observatory plus a signal from the 
Global Positioning System of satellites (GPS; [|13]). The TOA of this mean pulse is then 
defined as the arrival time of some fiducal point on the profile. Since the mean profile 
has a stable form at any given observing frequency (§ |2.2| ), the TOA can be accurately 
determined by a simple cross-correlation of the observed profile with a high signal-to-noise 
"template" profile obtained from the addition of many observations of the pulse profile 
at the particular observing frequency. 

Success in pulsar timing hinges on how precisely the fiducial point can be determined. 
This is largely dependent on the signal-to-noise ratio (SNR) of the mean pulse profile. 
The uncertainty in a TOA measurement e-roA is given roughly by the pulse width divided 
by the SNR. Using Equation ^, we can express this as a fraction of the pulse period: 

In this expression Tree and Tgky are the receiver and sky noise temperatures, G is the gain 
of the antenna, Au is the observing bandwidth, tint is the integration time, W is the pulse 
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width and P is the pulse period (we assume W ^ P). Optimum results are thus obtained 
for observations of short period pulsars with large flux densities and narrow duty cycles 
iW/P) using large telescopes with low-noise receivers and large observing bandwidths. 

One of the main problems of employing large bandwidths is pulse dispersion. As 
discussed in § |2.3| , the velocity of the pulsed radiation through the ionised interstellar 
medium is frequency-dependent: pulses emitted at higher radio frequencies travel faster 
and arrive earlier than those emitted at lower frequencies. This process has the effect of 
"stretching" the pulse across a finite receiver bandwidth, reducing the apparent signal- 
to- noise ratio and therefore increasing exoA- For most normal pulsars, this process can 



largely be compensated for by the incoherent de-dispersion process outlined in $^1 . 

To exploit the precision offered by millisecond pulsars, a more precise method of dis- 
persion removal is required. Technical difficulties in building devices with very narrow 
channel bandwidths require another dispersion removal technique. In the process of co- 
herent de-dispersion ||108|| the incoming signals are de-dispersed over the whole bandwidth 



using a filter which has the inverse transfer function to that of the interstellar medium. 
The signal processing can be done either on-line using finite impulse response filter devices 
|33| or off-line in software [p26| , |229|] . The on-line approach allows for large bandwidths to 



be employed and real-time viewing of the data. Off-line reduction, while slow and compu- 
tationally expensive, allows for more flexible data reduction schemes as well as periodicity 
searches to be carried out. 

The maximum time resolution obtainable via coherent dedispersion is the inverse of 
the receiver bandwidth. For bandwidths of 10 MHz, this technique makes it possible to 
resolve features on time-scales as short as 100 ns. This corresponds to probing regions in 
the neutron star magnetosphere as small as 30 m! 

4.2 The timing model 

Ideally, in order to model the rotational behaviour of the neutron star, we require TOAs 
measured by an inertial observer. An observatory located on Earth experiences acceler- 
ations with respect to the neutron star due to the Earth's rotation and orbital motion 
around the Sun and is therefore not in an inertial frame. To a very good approximation, 
the centre-of-mass of the solar system, the solar system barycentre, can be regarded as 



an inertial frame. It is standard practice ||120|| to transform the observed TOAs to this 



frame using a planetary ephemeris such as the JPL DE200 [231 1. The transformation is 



summarised as the difference between barycentric (T) and observed (t) TOAs: 

T - t = =:= + '-' + Atrci - AtDM. (8) 

c 2ca 

Here r is the position of the Earth with respect to the barycentre, § is a unit vector 
in the direction towards the pulsar at a distance d, and c is the speed of light. The 
first term on the right hand side of this expression is the light travel time from the 
Earth to the solar system barycentre. For all but the nearest pulsars, the incoming 
pulses can be approximated by plane wavefronts. The second term, which represents 
the delay due to spherical wavefronts and which yields the trigonometric parallax and 
hence d, is presently only measurable for four nearby millisecond pulsars ||134| , |65| , |215 . 



The term Atj-ei represents the Einstein and Shapiro corrections due to general relativistic 
effects within the solar system [^. Since measurements are often carried out at different 
observing frequencies with different dispersive delays, the TOAs are generally referred to 
the equivalent time that would be observed at infinite frequency. This transformation 
corresponds to the term A^dm and may be calculated from Equation ^ 
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Following the accumulation of about ten to twenty barycentric TOAs from observa- 
tions spaced over at least several months, a surprisingly simple model can be applied to 
the TOAs and optimised so that it is sufficient to account for the arrival time of any 
pulse emitted during the time span of the observations and predict the arrival times of 
subsequent pulses. The model is based on a Taylor expansion of the angular rotational 
frequency fl = 27i/P about a model value flo at some reference epoch to- The model pulse 
phase as a function of barycentric time is thus given by: 

0(T) = 0, + (T - To)^]o + ^(T - Tofdo + ■ • • , (9) 

where (po is the pulse phase at Based on this simple model, and using initial estimates 
of the position, dispersion measure and pulse period, a "timing residual" is calculated for 
each TOA as the difference between the observed and predicted pulse phases. 

A set of timing residuals for the nearby pulsar B1133+16 spanning almost 10 years 
is shown for illustrative purposes in Fig. [1^. Ideally, the residuals should have a zero 
mean and be free from any systematic trends (Fig. |IU|a). Inevitably, however, due to 
our a-priori ignorance of the rotational parameters, the model needs to be refined in a 
bootstrap fashion. Early sets of residuals will exhibit a number of trends indicating a 
systematic error in one or more of the model parameters, or a parameter not initially 
incorporated into the model. From Equation ^ an error in the assumed Qo results in a 
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Figure 19: Timing model residuals versus date for PSR B1133+16. Case (a) shows the 
residuals obtained from the best fitting model which includes period, the period derivative, 
position and proper motion. Case (b) is the result of setting the period derivative term 
to zero in this model. Case (c) shows the effect of a 1 arcmin error in the assumed 
declination. Case (d) shows the residuals obtained assuming zero proper motion. The 
lines in (b)-(d) show the expected behaviour in the residuals for each effect (see text). 



linear slope with time. A parabolic trend results from an error in fio (Fig. [19|b). Additional 
effects will arise if the assumed position of the pulsar (the unit vector s in equation ^) used 
in the barycentric time calculation is incorrect. A position error of just one arcsecond 
results in an annual sinusoid (Fig. P^Pjc) with a peak-to-peak amplitude of about 5 ms for a 
pulsar on the ecliptic; this is easily measurable for typical TOA uncertainties of order one 
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milliperiod or better. A proper motion produces an annual sinusoid of linearly increasing 
magnitude (Fig. p!9|d). 

After a number of iterations, and with the benefit of a modicum of experience, it is 
possible to identify and account for each of these various effects to produce a "timing 
solution" which is phase coherent over the whole data span. The resulting model parame- 
ters provide spin and astrometric information about the neutron star to a precision which 
improves as the length of the data span increases. Timing observations of the original 
millisecond pulsar, B1937-I-21, spanning almost 9 yr (exactly 165,711,423,279 rotations!) 
measure a period of 1.5578064688197945 ± 0.0000000000000004 ms [|133, |l3l defined at 
midnight UT on December 5 1988! Astrometric measurements based on these data are 
no less impressive, with position errors of ~ 20 /iarcsec being presently possible. 



4.3 Timing stability 

Ideally, after correctly applying a timing model, we would expect a set of uncorrelated 
timing residuals scattered in a Gaussian fashion about a zero mean with an rms consistent 
with the measurement uncertainties. This is not always the case; the residuals of many 
pulsars exhibit a quasi-periodic wandering with time. A number of examples are shown 
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Figure 20: Examples of timing residuals for a number of normal pulsars. Note the varying 
scale on the ordinate axis, the pulsars being ranked in increasing order of timing ^'activity". 



in Fig. These are taken from the Jodrell Bank timing program |p23| ]. Such "timing 
noise" is most prominent in the youngest of the normal pulsars [ 178 , 78|1 and virtually 
absent in the much older millisecond pulsars | |134| . While the physical processes of this 
phenomena are not well understood, it seems likely that it may be connected to superfluid 
processes and temperature changes in the interior of the neutron star |23] or processes in 



the magnetosphere |74, 73 



The relative dearth of timing noise for the older pulsars is a very important finding. 
It implies that, presently, the measurement precision depends primarily on the particular 
hardware constraints of the observing system. Consequently, a large effort in hardware 
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development is presently being made to improve the precision of these observations using, 
in particular, coherent dedispersion outlined in §[4.1|. Much of the pioneering work in 



this area has been made by Joseph Taylor and collaborators at Princeton University 
From high quality observations made using the Arecibo radio telescope spanning almost a 
decade [pl3| , pl4| , |134|| , the group has demonstrated that the timing stability of millisecond 
pulsars over such time-scales is comparable to terrestrial atomic clocks. 
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Figure 21: Fractional timing instabilities for PSRs B1855+09 and B1937+21 as a function 
of time. (After Kaspi, Taylor & Ryba 1994 /|7^/j. 



This phenomenal stability is demonstrated in Fig. P^. This figure shows a^, a param- 
eter closely resembling the Allan variance used by the clock community to estimate the 
stability of atomic clocks p38| , Atomic clocks are known to have ~ 5 x 10~^^ on 



time-scales of order 5 years. The timing stability of PSR B1937+21 seems to be limited 
by a power law component which produces a minimum in its az after ~ 2 yr. This is most 
likely a result of a small amount of intrinsic timing noise ||134|| . No such noise component 
is observed for PSR B1855+09. This demonstrates that the timing stability for PSR 
B1855-I-09 becomes competitive with the atomic clocks after about 3 yr. The absence of 
timing noise for B1855-I-09 is probably related to its characteristic age ~ 5 Gyr which is 
about a factor of 20 larger than B1937-t-21. Timing observations of millisecond pulsars 
are discussed further in the context of the pulsar timing array in 



5.2. 



4.4 Binary pulsars and Kepler's laws 



For binary pulsars, the simple timing model introduced in §4.2 needs to be extended 



to incorporate the additional radial acceleration of the pulsar as it orbits the common 
centre-of-mass of the binary system. Treating the binary orbit using Kepler's laws to refer 
the TOAs to the binary barycentre requires five additional model parameters: the orbital 
period (Pb), projected semi-major orbital axis (opsini, see below), orbital eccentricity (e), 
longitude of periastron (a;) and the epoch of periastron passage (To). This description, 
using five "Keplerian parameters" , is identical to that used for spectroscopic binary stars. 

For spectroscopic binaries the orbital velocity curve shows the radial component of the 
star's velocity as a function of time. The analogous plot for pulsars is the apparent pulse 
period against time. Two examples are given in Fig. |22|. Constraints on the mass of the 
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Figure 22: Orbital velocity curves for two binary pulsars. Left: PSR J1012+5307, a 
5.25-nis pulsar in a 14.5-hour circular orbit around a low-mass white dwarf companion 
liM m, mi- Right: PSR J1811-1736, a 104-ms pulsar in a highly eccentric 18.8-day 
orbit around a massive companion (probably another neutron star) 



orbiting companion can be placed by combining the projected semi-major axis Op sin i 
and the orbital period to obtain the mass function: 



Att"^ (Op sin i)^ 



G 



p2 



[mr sm I 



(mp + m^cY 



(10) 



where G is the universal gravitational constant. Assuming a pulsar mass mp of 1.35 Mq 
(see below), the mass of the orbiting companion mc can be estimated as a function of 
the (initially unknown) angle i between the orbital plane and the plane of the sky. The 
minimum companion mass mmin occurs when the orbit is assumed edge-on [i = 90°). For 
a random distribution of orbital inclination angles, the probability of observing a binary 
system at an angle less than some value io is p{< io) = 1 — cos(zo)- This implies that the 
chances of observing a binary system inclined at an angle < 26° is only 10%; evaluating 
the companion mass for this inclination angle mgo constrains the mass range between 
mmin and mgo at the 90% confidence level. 



4.5 Post-Keplerian parameters 

Although many of the presently known binary pulsar systems can be adequately de- 
scribed by a Keplerian orbit, there are several systems, including the original binary 
pulsar B1913+16, which exhibit relativistic effects that require an additional set of up to 
five "post-Keplerian" parameters. Within the framework of general relativity, these can 
be written 1551 as: 
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r = T0mc, (14) 
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In addition to the symbols defined above for Equation [Ty, M = rup + rric, x = OpSini/c, 
s = sini and Tq = GMq/c' ~ 4.925 yus. All masses are in solar units. 

Measurements of post-Keplerian parameters for PSR B1913+16 have been carried out 
by Taylor and a number of collaborators over the years with steadily improving precision. 
The first of these parameters to be measured was the advance of the longitude of periastron 
{(jj). This measurement is analogous to the perihelion advance of Mercury [|182 |. For PSR 
B1913+16 this amounts to about 4.2 degrees per year ||242|| , some 4.6 orders of magnitude 
larger than for Mercury. A measurement of uj alone yields the total mass for this system, 
M = 2.83 Mq assuming this advance is due to general relativity. 

Measurement of a second post-Keplerian parameter for B1913+16, 7, (gravitational 
redshift and transverse Doppler shifts in the orbit) permits an unambiguous determination 
of rrip, rric and i when combined with u and the five Keplerian parameters. The original 
measurements [MT] have since been substantially refined | 244 , 253] and the mass of the 
pulsar and its unseen companion have been determined to be 1.442±0.003 Mq and 1.386± 
0.003 Mq respectively. Such phenomenal precision is a testament to the timing stability 
of radio pulsars as clocks, and the diligence of Taylor and collaborators in carrying out 
these long-term measurements. Similar mass measurements now exist for the two other 
double neutron star binary systems discussed in § 3.4.1 : B1534+12 227 | and B2127+11C 
For a number of other systems, uj measurements allow interesting constraints to be 



placed on the component masses ||248| , |189| , |250| 



An important general relativistic prediction for eccentric double neutron star systems 
is the orbital decay due to the emission of gravitational radiation (Pf, in Equation |T3|). 
Taylor et al. |P41 , |244| , 245 1 were able to measure this for B1913+16 and found it to be in 
excellent agreement with the predicted value. The orbital decay, which corresponds to a 
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Figure 23: Orbital decay in the binary pulsar B1913+16 system demonstrated as an 
increasing orbital phase shift for periastron passages with time. The general relativistic 
prediction due entirely to the emission of gravitational radiation is shown by the parabola. 

shrinkage of about 3.2 mm per orbit, is seen most dramatically as the gradually increasing 
shift in orbital phase for periastron passages with respect to a non-decaying orbit shown 
in Fig. ^ This figure includes recent Arecibo data taken in 1998 and 1999 following 
the upgrade of the telescope in the mid 1990s. The observations of the orbital decay, 
now spanning a 25-year baseline, are in agreement with general relativity at the level 
of about 0.5% and provide the first (indirect) evidence for the existence of gravitational 
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radiation. Hulse and Taylor were awarded the Nobel prize in Physics in 1993 ||T], 118| , |239|| 
in recognition of their discovery of this remarkable laboratory for testing general relativity. 

For those binary systems which are oriented nearly edge-on to the line-of-sight, a 
significant delay is expected for orbital phases around superior conjunction where the 
pulsar radiation is bent in the gravitational potential well of the companion star. The 
so-called "range" and "shape" of the Shapiro delay effect are parameterized by the last 
two post-Keplerian parameters r and s = sin? that were introduced in Equations |1^ and 
151 . This effect, analogous to the solar system Shapiro delay, has so far been measured for 
two neutron star-white dwarf binary systems: B1855-I-09 and J1713+0747 pi3|, |134|, |65|] 



and for the double neutron star binaries B1534+12 227 and B1913+16 245 



4.6 Geodetic precession 



Shortly after the discovery of PSR B1913+16 it was realized that, if the spin axis of 
the visible pulsar was misaligned with the angular momentum axis of the binary system, 
the perturbing effect of the companion on the space-time around the radio pulsar would 



cause it to precess around the angular momentum axis 0. Within the framework of 



general relativity, the rate of precession ilp was shown |^ to be 



P^/' {mp + m,)4/3 (1 _ e2) 



(16) 



where we assume the same notation used for the discussion in § [4.4| and § [4.5[ Inserting the 
parameters of PSR B1913+16 yields Qp = 1.21 deg yr~^. The period of the precession is 
297.5 yr. The observational consequence of geodetic precession is a secular change in the 
pulse profile as the line-of-sight cut through the emission beam changes (recall Fig. 





7975 



Figure 24: Geodetic precession in the binary pulsar 31913+16 system: (a) changes in the 
observed pulse shapes between 1981-1988 seen as a decrease in amplitude between the 
left and right components in the pulse profile \24d[ , ^66^ ; (b) relative heights of the two 
amplitudes plotted between 1975-1998 JJ36|/ ; (c) component separation change ^136^] . 



Early qualitative evidence for profile evolution due to this effect p41|| was substantiated 
with long-term Arecibo measurements of component changes by Weisberg et al. ||266| | . 
Further changes were seen by Kramer with new Effelsberg data acquired in the 1990s 
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136 1 . In addition to relative amplitude variations, the expected changes in component 
separation for a hollow-cone beam model were also seen in the Effelsberg data. These 



observations are summarized in Fig. 

In addition to the above results, there is now evidence for geodetic precession in the 
other classic neutron star binary, PSR B1534+12 p9| , |23Q|| . Although geodetic precession 
in binary pulsars is another successful test of general relativity (albeit at a lower precision 
than e.g. orbital decay measurements), what is perhaps more interesting are the various 
consequences it has. Geodetic precession only occurs when the spin and orbital axes are 
misaligned This is most likely to occur if the neutron star received an impulsive 

"kick" velocity at birth (§|A|). Wex et al. [^E^ have investigated the B1913+16 obser- 
vations and find that the kick magnitude was at least 250 km s~^ and was directed almost 
perpendicular to the spin axis of the neutron star progenitor. This places stringent con- 
straints on any kick mechanism. Detailed monitoring of the pulse profile and polarization 
properties now underway [p64| , |137|| will allow the first map of the emission beam of a 
neutron star to be made. This has important implications for the various beaming models 
described in 



3.2.3 



There are already indications that the beam is circular ||137|| . 
The current results predict that B1913+16 will completely precess out of the line 
of sight by around 2025 and re-appear some 240 years later |136|. Although we shall 



lose a most treasured pulsar, we can take comfort from the fact that other pulsars will 
precess into our field of view. Perhaps one example is the newly-discovered relativistic 
binary J1141— 6545 [ |133|| discussed in § p.6.2| and § |3.4.2| . This relatively bright object was 
apparently missed by two previous searches during the early 1990s ||125| , |176| , |168|| . 



4.7 Going further 

This chapter has outlined past and present progress in a number of areas related to pulsar 
timing. For further details on the technical details and prospects of pulsar timing, the 
interested reader is referred to a number of excellent review articles ||270| , |35i 



Two freely available software packages which are routinely used for time-of-arrival analyses 
by the pulsar community are available viz: TEMPO p45|, |] and TIMAPR 0, [l^. 
These packages are based on more detailed versions of the timing model outlined in 
An up-to-date list summarising the various timing programmes is kept by Don 
An audio file and slides from a lecture on pulsar timing presented by Backer 



4.2 



Backer | 

at the centennial meeting of the American Physical society is also available on-line 



Another relevant lecture from that meeting is Will's presentation on tests of Einstein's 
relativity which includes an excellent overview of Taylor and Weisberg's measurements of 
PSR B1913-I-16. Kramer ||137|| has written a lucid review article discussing measurements 
of geodetic precession in binary pulsars and their implications. 

The remarkable precision of these measurements, particularly for millisecond pulsars, 
allows the detection of radial accelerations on the pulsar induced by orbiting bodies smaller 
than the Earth. Alex Wolszczan detected one such "pulsar planetary system" in 1990 fol- 
lowing the discovery of a 6.2-ms pulsar B1257-I-12. In this case, the pulsar is orbited by at 
least three Earth-mass bodies ||274|, E 



^ ^ 273|| . Subsequent measurements of B1257-I-12 were 

even able to measure resonance interactions between two of the planets [ f272|| , confirming 
the nature of the system beyond all doubt. Long-term timing measurements of the 11-ms 
pulsar B1620-26 in the globular cluster M4 indicate that it may also have a planetary 
companion p49| , |3^ , p47|| . For detailed reviews of these systems, and their implications 



for planetary formation scenarios, the interested reader is referred to [198 . 
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5 Pulsars as Gravitational Wave Detectors 



Many cosmological models predict that the Universe is presently filled with a stochastic 
gravitational wave background (GWB) produced during the big bang era [|195|| . The idea 
to use pulsars as natural detectors of gravitational waves was first explored independently 
by both Sazhin and Detweiler in the late 1970s ||216| , [86|] . The basic concept is to treat 
the solar system barycentre and a distant pulsar as opposite ends of an imaginary arm in 
space. The pulsar acts as the reference clock at one end of the arm sending out regular 
signals which are monitored by an observer on the Earth over some time-scale T. The 
effect of a passing gravitational wave would be to cause a change in the observed rotational 
frequency by an amount proportional to the amplitude of the wave. For regular monitoring 
observations of a pulsar with typical TOA uncertainties of exoA, this "detector" would 
be sensitive to waves with dimensionless amplitudes e^ox/T and frequencies as low as 



1/T [mn 



is reviewed in S5.1 



This method, which already yields interesting upper limits on the GWB, 
The idea of more sensitive detector based on an array of pulsar clocks 



distributed over the sky is discussed in §5.2 



5.1 Limits from individual pulsars 

In the ideal case, the change in the observed frequency caused by the GWB should be 
detectable in the set of timing residuals after the application of an appropriate model 
for the rotational, astrometric and, where necessary, binary parameters of the pulsar. 
As discussed in §§, all other effects being negligible, the rms scatter of these residuals 
a would be due to the measurement uncertainties and intrinsic timing noise from the 



neutron star. Detweiler p6[ showed that a GWB with a flat energy spectrum in the 
frequency band / ± //2 would result in an additional contribution to the timing residuals 
ag. The corresponding wave energy density pg (for /T ^ 1) is 

243 7r3/'(T2 ^ ^ 

= 208 G ■ ^'^^ 

An upper limit to Pg can be obtained from a set of timing residuals by assuming the rms 
scatter is entirely due to this effect (cr = (Xg). These limits are commonly expressed as a 
fraction of pc the energy density required to close the Universe: 

3 

p, = ~ 2 X 10-29 g cm-3, (18) 

where the Hubble constant Ho = 100 h km s~^ Mpc. 

Romani & Taylor |pi2|| applied this technique to a set of TOAs for PSR B1237+12 



obtained from regular observations over a period of 11 years as part of the JPL pulsar 



timing programme |91|. This pulsar was chosen on the basis of its relatively low level of 
timing activity by comparison with the youngest pulsars, whose residuals are ultimately 
plagued by timing noise (§ [4.3|) . By ascribing the rms scatter in the residuals {a = 240 
ms) to the GWB, Romani & Taylor placed a limit of Pg/pc ~ 4 x 10~^h~'^ for a centre 
frequency / = 7 x 10"^ Hz. 

This limit, already well below the energy density required to close the Universe, was 
further reduced following the long-term timing measurements of millisecond pulsars at 
Arecibo by Taylor and collaborators (§ [4.3|) . In the intervening period, more elaborate 
techniques had been devised ^34|| to look for the likely signature of a GWB in the 



frequency spectrum of the timing residuals and to address the possibility of "fltting out' 
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the signal in the TOAs. Following |H8| it is convenient to define fig, the energy density 



of the GWB per logarithmic frequency interval relative to pc- With this definition, the 



power spectrum of the GWB, V{f), can be written [|117| , |56 



as 



G Pg _ Hi Vtg 

37r3 /4 ~ 8774 /5 



1.34X 10^^]3/^Vy",-l /WsVr, 



(19) 



where /yr-i is frequency in cycles per year. The timing residuals for B1937+21 shown in 
Fig. ^ are clearly non- white and, as we saw in §0, limit its timing stability for periods 
^ 2 yr. The residuals for PSR B 1855+09 clearly show no systematic trends and are in 
fact consistent with the measurement uncertainties alone. Based on these data, and using 
a rigorous statistical analysis, Thorsett & Dewey ||251|| place a 95% confidence upper limit 
of flgh? < 10~^ for / = 4.4 x 10~^ Hz. This limit is difficult to reconcile with most cosmic 



string models for galaxv formation lEPI. 1251 1. 
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Figure 25: Timing residuals for PSRs B1855+09 (left panel) 31937+21 and (right panel) 
obtained from almost a decade of timing at Areciho. (Kaspi, Taylor & Ryha 1994 JJ34[ /j. 

For those pulsars in binary systems, an additional clock for measuring the effects of 
gravitational waves is the orbital period. In this case, the range of frequencies is not 
limited by the time span of the observations, allowing the detection of waves with periods 



as large as the light travel time to the binary system [Q. The most stringent results 
presently available are based on B1855+09 limit flgh"^ < 2.7 x 10^^ in the frequency range 
10~^^ < / < 4.4 X 10~^ Hz. Kopeikin ||135|| has recently presented this limit and discusses 
the methods in detail. 



5.2 A pulsar timing array 

The idea of using timing data for a number of pulsars distributed on the sky to detect 



gravitational waves was first proposed by Hellings & Downs [p.l4|| . Such a "timing array" 
of pulsars would have the advantage over a single arm in that, through a cross- correlation 
analysis of the residuals for pairs of pulsars distributed over the sky, it should be possible 
to separate the timing noise of each pulsar from the signature of the GWB, which would 
be common to all pulsars in the array. To quantify this, consider the fractional frequency 
shift observed for the i^^ pulsar in the array: 

^ = a,^(t)+M(t). (20) 

In this expression a, is a geometric factor dependent on the line-of-sight direction to 
the pulsar and the propagation and polarisation vectors of the gravitational wave of 
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dimensionless amplitude A. The timing noise intrinsic to the pulsar is characterised by 
the function A/^. The result of a cross-correlation between pulsars i and j is then 



> +ai < AAfj > +aj < AAfi > + < UiUj >, 



f211 



where the bracketed terms indicate cross- correlations. Since the wave function and the 
noise contributions from the two pulsars are independent quantities, the cross correlation 
tends to aiaj < > as the number of residuals becomes large. Summing the cross- 
correlation functions over a large number of pulsar pairs provides additional information 
on this term as a function of the angle on the sky ||113|| . This allows the separation of the 



effects of terrestrial clock and solar system ephemeris errors from the GWB |97 



Applying the timing array concept to the present database of long-term timing ob- 
servations of millisecond pulsars does not improve on the limits on the GWB discussed 
above. The sky distribution of these pulsars, seen in the left panel of Fig. |26|, shows 
that their angular separation is rather low. To achieve optimum sensitivity it is desirable 
to have an array consisting of pulsar clocks distributed isotropically over the whole sky. 
The flood of recent discoveries of nearby binary and millisecond pulsars has resulted in 
essentially such a distribution, shown in the right panel of Fig. |2^. 



1990 



2001 





Figure 26: Hammer- AitoS projections showing the known Galactic disk miUisecond pulsar 
population in 1990 and 2001. The impact of the new discoveries is seen by comparing the 
sample circa 1990, where all the known sources had been discovered at Arecibo, with the 
present sample where the sources are much more uniformly distributed on the sky 

A number of long-term timing projects are now underway to monitor these millisec- 
ond pulsars with a goal of detecting low-frequency gravitational radiation. At Arecibo, 
regular timing of a dozen or more millisecond pulsars has been carried out following the 
completion of the upgrade to the telescope in 1997. A summary of these observations is 



shown in Fig. ^7} The rms timing residuals for several of the pulsars are now approach- 
ing the 100 ns level. This degree of precision demands a high level of commitment to 
investigate possible causes of systematic errors in the signal path through the telescope. 
Combining datasets from several observatories is also challenging. The Berkeley pulsar 
group lead by Don Backer ^ are among the most active observers in this area. Backer and 
collaborators have now installed identical sets of datataking equipment at a number of 
radio telescopes around the world in an attempt to ensure a homogeneous set of residuals. 
Continued timing of these millisecond pulsars in the coming years should greatly improve 
the sensitivity and will perhaps allow the detection of gravitational waves, as opposed to 
upper limits, in the not-too-distant future. 
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Figure 27: Arecibo timing residuals of millisecond pulsars monitored by groups at Berkeley 
and Princeton on a regular basis as part of a long-term project to detect low- frequency 
gravitational waves. The level of precision being achieved is well below the vertical bars 
surrounding each set of residuals (±5fis). 



5.3 Going further 



Further discussions on the reahties of using pulsars as gravity wave detectors can be found 
in two excellent review articles by Romani | |210| | and Backer . 
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6 Summary and Future Prospects 



The main aim of this article was to review some of the many astrophysical apphcations 
provided by the present sample of binary and millisecond radio pulsars. The topics cov- 
ered here, along with the bibliography and associated tables of observational parameters, 
should be useful to those wishing to delve deeper into the vast body of literature that 
exists. We now briefly recap on the main issues. 

Through an understanding of the Galactic population of radio pulsars summarised in 
§^ it is possible to predict the detection statistics of terrestrial gravitational wave detectors 
to nearby rapidly spinning neutron stars (§ |3.3| ), as well as coalescing relativistic binaries at 
cosmic distances (§ |3.4|) . Continued improvements in gravitational wave detector sensitiv- 
ities should result in a number of interesting developments and contributions in this area. 
These developments and contributions might include the detection of presently known ra- 
dio pulsars, as well as a population of coalescing binary systems which have not yet been 
detected as radio pulsars. The phenomenal timing stability of radio pulsars leads natu- 
rally to a large number of applications, including their use as laboratories for relativistic 
gravity (§ [4.5| ) and as natural detectors of gravitational radiation (§^. Long-term timing 



experiments of the present sample of millisecond and binary pulsars currently underway 
appear to have tremendous potential in these areas and perhaps detect the gravitational 
wave background (if it exists) within the next decade. 

These applications will benefit greatly from the continued discovery of new systems 
by the present generation of radio pulsar searches which continue to probe new areas of 
parameter space. Based on the results presented in $3. 3 , it is clear that we are aware of 
only about 1% of the total active pulsar population in our Galaxy. It is therefore likely 
that we have not seen all of the pulsar zoo. More sensitive surveys are being planned 
both in the short term (a multibeam system on the Arecibo telescope [|18|) and in the 
longer term (the Square Kilometer Array |jl9|)- These should provide a far more complete 



census of the Galactic pulsar population. Possible discoveries in the future include: 

• A dual-line binary pulsar, i.e. a double neutron star system in which both compo- 
nents are observable as radio pulsars. The additional clock in such a binary system 
would be most valuable in further tests of strong-field gravity. 

• A radio pulsar with a black-hole companion would undoubtably also be a fantastic 
laboratory for studying gravity in the strong-field regime. 

• A sub-millisecond pulsar. The original millisecond pulsar, B1937-I-21, rotating at 
642 Hz is still the most rapidly rotating neutron star known. Do kHz neutron stars 
exist? Searches now have sensitivity to such objects and a discovery of even 



one would constrain the equation of state of matter at high densities. 

A binary system in which the neutron star is in the process of transforming from 
an X-ray-emitting neutron star to a millisecond radio pulsar. 



Spurred on by recent discoveries p6|, |1 70|, |206|, E^ , a number of high-sensitivity searches 



for pulsars in globular clusters are being conducted. These have tremendous potential for 
discovering new and exotic binary systems like a millisecond pulsar-black hole binary. 
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A Tables of Binary and Millisecond Pulsars 



Name 



P log(r,) log(B) 



Vt 



Ref. 







ms 






kpc 


km/s 


J00304 


-0451 


4.865 


9.9 


8.3 


0.23 


<65 


J0711- 


6830 


5.491 


10.4 


8.2 


1.04 


139 


J1024- 


0719 


5.162 


9.7 


8.5 


0.35 


45 


J1730- 


2304 


8.123 


9.9 


8.6 


0.51 


53 


J1744- 


-1134 


4.075 


9.9 


8.3 


0.17 


20 


B1937- 


h21 


1.558 


8.4 


8.6 


9.65 


22 


J2124- 


-3358 


4.931 


9.9 


8.4 


0.25 


67 


J22354 


-1506 


59.767 


10.0 


9.4 


1.15 


98 


J23224 


-2057 


4.808 


10.2 


8.2 


0.78 


89 




Table 2: Parameters for the 9 isolated miUisecond pulsars currently known in the Galactic 
disk. Listed are the spin period P, the base-10 logarithms of the characteric age Tc and 
surface magnetic field strength B ("§ |2.4. the distance d derived from the Taylor <fe Cordes 
electron density model ^240^ or independently (when available) and the transverse speed 
Vt inferred from d and a proper motion measurement (when available). Key publications 
for each pulsar are referenced to the bibliography. 
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Name 


P log(re) log(B) d 




n 


X 


e TO2 


Ref. 




ms kpc 


km/s 


days 


s 


Mo 





J1141-6545 

B1259-63 

J1518+4904 

B1534+12 

J1740-3052 

J1811-1736 

B1820-11 
B1913+16 

B2303+46 



393.898 
47.762 
40.935 

37.904 
570.309 
104.182 

279.828 
59.030 
1066.371 



6.2 
5.5 
10.3 

8.4 
5.5 
8.9 

6.5 
8.0 
7.5 



12.1 
11.5 
9.0 

10.0 
12.6 
10.1 

11.8 
10.4 
11.9 



3.20 
4.60 
0.70 

0.68 
10.8 
5.94 

6.26 
7.13 
4.35 



? 
? 
27 

80 

? 
? 

? 

100 

? 



0.20 
1236.72 
8.63 

0.42 
231.03 
18.77 

357.76 
0.32 
12.34 



1.86 
1296.58 
20.04 

3.73 
756.91 
34.78 

200.67 
2.34 
32.69 



0.17 
0.87 
0.25 

0.27 
0.58 
0.83 

0.79 
0.62 
0.66 



1.0 
10.0 
1.3 

1.3 
16 
0.7 

0.7 
1.4 
1.2 



125 
189 



13S 



12f I 
19C| 



271 



22/ 



|163| 



171 

244 
159 



202 1 



245 1 



25^ 



Table 3: Parameters for the 9 high-eccentricity (e > 0.15) binary pulsars currently known 
in the Galactic disk. Listed are the spin period P, the base-10 logarithms of the characteric 
age Tc and surface magnetic field strength B (^§ |2.4.i[ ), the distance d derived from the 
Taylor & Cordes electron density model ^240^ or independently (when available), the 
transverse speed Vt inferred from d and a proper motion measurement (when available), 
the binary period Pb, the projected semi-major axis of the orbit x in units of light seconds, 
the orbital eccentricity e and the companion mass m2 evaluated from the mass function 
assuming a pulsar mass of 1.4 Mq and an inclination angle of 60 degrees f§ [4.4 ) or (when 
known) from independent measurements. Key publications for each pulsar are referenced 
to the bibliography. 
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Name 


P log(T,) log(B) 


d 




Pb 


X 


e 1112 


Ref. 




ms 


kpc 


km/s 


days 


s 


Mo 





J0034-0534 


1 


877 


9 


9 


7 


9 





98 


71 


1 


59 


1.44 


< 0.00002 


0.1 


J0218+4232 


2 


323 


8 


7 


8 


6 


5 


85 


? 


2 


03 


1.98 


< 0.00002 


0.2 


J0437-4715 


5 


757 


9 


7 


8 


5 





18 


121 


5 


74 


3.37 


0.000019 


0.1 


J0613-0200 


3 


062 


9 


7 


8 


2 


2 


19 


77 


1 


20 


1.09 


0.000007 


0.1 


J0621+1002 


28 


854 


10 


1 


9 





1 


88 


? 


8 


31 


12.03 


0.0025 


0.5 


130655+04 


195 


671 


9 


7 


10 


1 





A O 

48 


32 


1 


03 


4.13 


0.000008 


A 'T 

0.7 


J0751+1807 


3 


479 


9 


8 


8 


2 


2 


02 


? 





26 


0.40 


> 0.0 


0.1 


B0820+02 


864 


873 


8 


1 


11 


5 


1 


43 


35 


1232 


47 


162.15 


0.011868 


0.2 


J1012+5307 


5 


256 


9 


8 


8 


4 





52 


102 





60 


0.58 


< 0.0000008 


0.1 


J1022+1001 


16 


453 


9 


8 


8 


9 





60 


> 50 


7 


81 


16.77 


0.000098 


0.7 


J1045-4509 


7 


474 


10 





8 


5 


3 


25 


52 


4 


08 


3.02 


0.000024 


0.2 


J1157-5114 


43 


589 


9 


7 


9 


4 


1 


88 


? 


3 


51 


14.29 


0.00040 


1.2 


J1232-6501 


88 


282 


9 


2 


9 


9 


10 


00 


? 


1 


86 


1.61 


0.00011 


0.1 


B1257+12 


6 


219 


9 


5 


8 


6 





62 


281 






planetary system 




J1435-6100 


9 


348 


9 


8 


8 


7 


3 


25 


? 


1 


35 


6.18 


0.00001 


0.9 


J 1454— 5845 


45 


249 


9 





9 


o 

8 


3 


32 


o 

f 


12 


42 


26.52 


0.0019 


A A 

0.9 


J1455-3330 


7 


987 


9 


9 


8 


5 





74 


100 


76 


17 


32.36 


0.000170 


0.3 


J1603-7202 


14 


842 


10 


2 


8 


7 


1 


64 


27 


6 


31 


6.88 


> 0.0 


0.3 


J1640+2224 


3 


163 


10 


2 


8 





1 


18 


76 


175 


46 


55.33 


0.0008 


0.3 


J1643-1224 


4 


622 


9 


7 


8 


4 


4 


86 


159 


147 


02 


25.07 


0.000506 


0.1 


J1713+0747 


4 


570 


10 





8 


3 





89 


27 


67 


83 


32.34 


0.000075 


0.3 


J1757-5322 


8 


870 


9 


7 


8 


7 


1 


36 


? 





45 


2.09 


0.000004 


0.6 


B1800-27 


334 


415 


8 


5 


10 


9 


3 


62 


? 


406 


78 


58.94 


0.000507 


0.1 


J1804-2717 


9 


343 


9 


5 


8 


8 


1 


17 


? 


11 


13 


7.28 


0.000035 


0.2 


J1810-2005 


32 


822 


9 


6 


9 


3 


4 


04 


? 


15 


01 


11.98 


0.000025 


0.3 


rslool— UU 


con 


954 


o 



o 
<5 


1 A 
10 


9 


o 
Z 


63 




1 


81 


0.72 


> 0.0 


U.l 


B1855+09 


5 


362 


9 


7 


8 


5 


1 


00 


29.2 


12 


33 


9.23 


0.000022 


0.2 


J1904+0412 


71 


095 


10 


1 


9 


4 


4 


01 


? 


14 


93 


9.63 


0.00022 


0.2 


J1911-1114 


3 


626 














1 


59 


183 


2 


72 


1.76 


> 0.0 


0.1 


B1953+29 


6 


133 


9 


5 


8 


6 


5 


39 


98 


117 


35 


31.41 


0.00033 


0.2 


B1957+20 


1 


607 


9 


4 


8 


1 


1 


53 


190 





38 


0.09 


> 0.0 


0.02 


J2019+2425 


3 


935 


10 


4 


8 








91 


83 


76 


51 


38.77 


0.000111 


0.3 


J2033+1734 


5 


949 


9 


9 


8 


4 


1 


38 


? 


56 


31 


20.16 


0.00013 


0.2 


J2051-0827 


4 


509 


9 


7 


8 


4 


1 


28 


14 





10 


0.05 


> 0.0 


0.03 


J2129-5721 


3 


726 


9 


5 


8 


4 


2 


55 


56 


6 


63 


3.50 


> 0.0 


0.1 


J2145-0750 


16 


052 


10 


3 


8 


6 





50 


38 


6 


84 


10.16 


0.000019 


0.4 


J2229+2643 


2 


978 


10 


4 


7 


9 


1 


43 


113 


93 


02 


18.91 


0.00026 


0.1 


J2317+1439 


3 


445 


10 


6 


7 


9 


1 


89 


68 


2 


46 


2.31 


> 0.0 


0.2 



[152, 252| 



|185 
|124| 



[128, 260| 



177, 260 



[|187 
19 



[1561 



[152 



[274, 



[152 



[|151 



[|152 



[2731 



[1251 



[1511 



[221 



[ |151 , 252| 

^3 iz^i 



[101 



[192 



[232 



[151 



[208| 



[70 



[2731 



Table 4: Parameters for 38 low- eccentricity binary pulsars currently known in the Galactic 
disk. Listed are the spin period P, the base-10 logarithms of the characteric age Tc and 
surface magnetic field strength B ("§ 12.4. ip , the distance d derived from the Taylor & Cordes 

or independently (when available), the transverse speed Vt 



electron density model / |24( 
inferred from d and a proper motion measurement (when available), the binary period 
Pb, the projected semi-major axis of the orbit x in units of light seconds, the orbital 
eccentricity e (when measured) and the companion mass m2 evaluated from the mass 
function assuming a pulsar mass of 1.4 Mq and an inclination angle of 60 degrees ("§ 14. ■^ j 
or (when known) from independent measurements. Key publications for each pulsar are 
referenced to the bibliography. 
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PbR 


F Cluster 


7 

a 


^^6 




e 


m2 


Ret. 




(ms) 


(kpc) 


(days) 






Mo 





J0023-7204C 


5 


757 


47 Tuc 


4 


5 










J0024-7204D 


5 


358 


47 Tuc 


4 


5 










J0024-7205E 


3 


536 


47 Tuc 


4 


5 


2.26 


1.98 


0.0003 


0.2 


J0024-7204F 


2 


624 


47 Tuc 


4 


5 










J0024-7204G 


4 


040 


47 Tuc 


4 


5 










J0024-7204H 


3 


210 


47 Tuc 


4 


5 


2.36 


2.15 


0.07 


0.2 


J0024-7204I 


3 


485 


47 Tuc 


4 


5 


0.23 


0.04 


< 0.001 


0.01 


J0023-7203J 


2 


101 


47 Tuc 


4 


5 


0.12 


0.04 


< 0.0002 


0.02 


J0024-7204L 


4 


346 


47 Tuc 


4 


5 










J0023-7205M 


3 


677 


47 Tuc 


4 


5 










J0024-7204N 


3 


054 


47 Tuc 


4 


5 










J0024-7204O 


2 


643 


47 Tuc 


4 


5 


0.14 


0.05 


< 0.004 


0.02 


J0024-72P 


3 


643 


47 Tuc 


4 


5 


0.14 


0.04 


> 0.0 


0.02 


J0024-7204Q 


4 


033 


47 Tuc 


4 


5 


1.19 


1.46 


0.00007 


0.2 


J0024-72R 


3 


480 


47 Tuc 


4 


5 


0.066 


0.033 


> 0.0 


0.03 


J0024-72S 


2 


830 


47 Tuc 


4 


5 


1.20 


0.77 


> 0.0 


0.09 


J0024-7204T 


7 


589 


47 Tuc 


4 


5 


1.13 


1.34 


0.0004 




J0024-7203U 


4 


343 


47 Tuc 


4 


5 


0.43 


0.53 


0.0002 




J0024-72V 


4 


810 


47 Tuc 


4 


5 


curre 


ntly no 


orbital solution 


J0024-72W 


2 


352 


47 Tuc 


4 


5 


0.11 


0.24 


> 0.0 


0.15 


B1310+18 


33 


163 


M53 


18 


9 


55.80 


84.20 


0.002 


0.3 


B1516+02A 


5 


554 


M5 


7 













B1516+02B 


7 


947 


M5 


7 





6.86 


3.04 


0.14 




B1620-26 


11 


076 


M4 


1 


8 


91.44 


64.81 


0.03 


0.3 


B1639+36A 


10 


378 


M13 


7 


7 










B1639+36B 


3 


528 


M13 


7 


7 


1.26 


1.39 


0.005 


0.2 


J1701-30 


5 


242 


NGC6266 


6 


7 


3.81 


3.48 


> 0.0 


0.2 


B1718-19 


1004 


037 


NGC6342 


7 





0.26 


0.35 


> 0.0 


0.1 


J1740-53 


3 


650 


NGC6397 


2 


2 


1.35 


1.66 


> 0.0 


0.2 


B1744-24A 


11 


563 


Terzan5 


7 


1 


0.08 


0.12 


> 0.0 


0.1 


J1748-2446C 


8 


436 


Terzan5 


7 


1 










B1745-20 


288 


603 


NGC6440 


5 


8 










B1802-07 


23 


101 


NGC6539 


3 


1 


2.62 


3.92 


0.21 


0.3 


J1807-24 


3 


059 


NGC6544 


2 


5 


0.071 


0.012 


> 0.0 


0.009 


B1820-30A 


5 


440 


NGC6624 


8 













B1820-30B 


378 


596 


NGC6624 
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B1821-24 


3 
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5 
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3 
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3 


9 


0.865 


1.27 


> 0.0 


0.19 


J1910+0004 


3 


619 
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1 


0.14 


0.04 


> 0.0 


0.02 
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110 
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56 
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30 


529 


M15 


10 
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0.9 
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4 
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10 













B2127+11E 


4 


651 
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10 













B2127+11F 
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027 
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10 
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37 
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10 
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743 
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10 
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pi] 

11 
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17C1 



16fl 



24i 
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1 275 



J 25 



1 204 



1 204 



181] 
185] 
24] 
24] 
24] 
24] 
24] 

m 

\24] 
\24] 



Table 5: Parameters for the 47 pulsars currently known in globular clusters. Listed are 
the spin period P, cluster name and the distance d to the cluster. For binary pulsars we 
also list the binary period Pb, the projected semi-major axis of the orbit x in units of light 
seconds, the orbital eccentricity e and the companion mass m2 evaluated from the mass 
function assuming a pulsar mass of 1.4 Mq and an inclination angle of 60 degrees ( ^.4 ) 
or (when known) from independent measurements. Key publications for each pulsar are 
referenced to the bibliography. 
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